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^ : Abstract 

qh| a GHz maser beam with Gaussian-type distribution passing through a homoge- 

^ [ nous static magnetic field can be used to detect gravitational waves (GWs) with the 

same frequency. The presence of GWs will perturb the electromagnetic (EM) fields, 
giving rise to perturbed photon fluxes (PPFs). After being reflected by a fractal 
i/-^ ' membrane, the perturbed photons suffer little decay and can be measured by a 

microwave receiver. This idea has been explored to certain extent as a method for 
very high frequency gravitational waves. In this paper, we examine and develop 
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Q ■ this method more extensively, and confront the possible detection with the pre- 
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dieted signal of relic gravitational waves (RGWs). A maser beam with high linear 
polarization is used to reduce the background photon fluxes (BPFs) in the detect- 
K> , ing direction as the main noise. As a key factor of applicability of this method, 

j^ ■ we give a preliminary estimation of the sensitivity of a sample detector limited by 

thermal noise using currently common technology. The minimal detectable ampli- 
tude of GWs is found to be h^i^ ~ 10~^°. Comparing with the known spectrum of 
the RGWs in the accelerating universe for j3 = —1.9, there is still roughly a gap of 
4 ~ 5 orders. However, possible improvements on the detector can further narrow 
down the gap and make it a feasible method to detect high frequency RGWs. 
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I. INTRODUCTION 

Although there has been some indirect evidence of GWs radiation from the binary pulsar 
B1913+16 [1], so far GWs have not been directly detected yet. Currently, a number of detectors 
have been running or under construction. These detectors use various methods including: (1) 
the conventional method of cryogenic resonant bar [2] aiming at a frequency around 10^ Hz; 
(2) the method of ground-based laser interferometers, such as LIGO [3], VIRGO [4], GEO [5], 
TAMA [6], and AIGO [7], applying for a frequency range 10 ~ 10'^ Hz, and the space-based 
laser interferometers LISA [8] under planning for a lower frequency range (10~^, 10~^) Hz; (3) 
the detections of cosmic microwave background radiation (CMB) polarization of "magnetic" 
type, or the temperature- "electric" cross-correlation [9] , which would also give direct evidence 
of GWs [10] for very low frequencies around (10~^^, 10"^^) Hz. Besides, there have also been 
attempts, based on various techniques, to detect GWs of very high frequencies from MHz to 
GHz, such as the waveguide detector to measure the change of polarizations of EM waves 
[11, 12], the two coupled microwave cavities to measure small harmonic displacements [13], 
and laser interferometers [14], etc. 

There is another method of detection for high frequency gravitational waves (HFGWs), 
which employs a a maser beam passing through a strong static magnetic field [15, 16, 17], and 
uses a microwave receiver in combination with a fractal membrane [18, 19]. The maser beam 
can be chosen to a free electron maser [20, 21] with a great output power ~ 2 kW, whose 
frequency ~ 4.5 GHz is the one that the fractal membrane operates effectively [22, 23, 24]. In 
the presence of GWs, the background EM fields will be perturbed, giving rise to additional 
photon fluxes in various directions. In the previous studies by Li et al [17, 19], an ordinary 
maser beam is used, in which case the BPFs always exist and tend to mix up with the PPFs, 
forming a kind of noise. Moreover, in order to assess the method as a potential way to detect 
GWs, the sensitivity of the detection predicted by this method has to be estimated, and a 
comparison with the predicted spectrum of RGWs [25]-[29] is still needed. In this paper, 
we improve the method by using a linearly polarized maser beam, so that the BPFs in the 
detecting direction can be suppressed effectively. 

Generally speaking, HFGWs in GHz band are not generated by usual astrophysical pro- 
cesses, such as explosions of asymmetric supernovas, rotations of binary stars around each 
other, coalescing and merging of binary neutron stars or black holes, and collapse of stars 
[30, 31, 32]. There could be a thermal background of gravitational waves, which consists of 
gravitons in thermal equilibrium [33, 34]. But, it will be examined that, if the inflationary 
expansion did occur in the early universe, the possible thermal GWs will be negligibly small. 
As is known, RGWs generated by the inflation have a spectrum stretching over a whole range 
of (10^^^, 10^^) Hz. In particular, it has a considerable amount of power around the very high 
frequency range ~ 4.5 GHz. Thus it can serve as the main target of detection using the maser 
beam. We shall estimate qualitatively the sensitivity of detection of a sample detector, and 



make a comparison with the known analytic spectrum of RGWs [27, 28, 29], which has not 
been made before. 

The organization of the paper is as follows. In Section II, we describe a setup of maser- 
beam GW detection. In Section III we study the PPFs density caused by the incident GWs 
propagating along various directions, and estimate the number of the perturbed photons per 
second passing through a receiving surface. Section IV is devoted to a preliminary analysis 
of the sensitivity of the detector limited by the thermal noise, and to the discussion of the 
feasibility of detecting the RGWs in the accelerating universe. In Section V, a summary is 
given. The Appendix gives a detailed calculation of the PPFs density generated by the incident 
GWs along the positive z-direction. 

II. THE SETUP OF THE DETECTOR 

The idea of the maser beam gravitational wave detector is based on the property that the 
maser beam in the presence of a homogeneous static magnetic field will be perturbed when 
GWs pass by [16, 19]. In particular, under the resonance condition that the frequency of GWs 
equals that of the maser beam {ug = z/g), additional PPFs will be generated and serve as a 
signal of GWs to be detected. As can be seen later, the magnitude of the PPFs is proportional 
to the amplitude of GWs and to the static magnetic field as well. Fig.l shows the geometric 
configuration of the setup for the detector. The maser beam of frequency z/g ~ 4.5 GHz travels 
along the positive z-direction and passes through the static magnetic field ~ 3 Tesla pointing 
to the positive y-direction. The PPF density (photons per unit area per unit time) in the 
x-direction, n'^\ after being totally reflected by a fractal membrane [22, 23, 24], will keep its 
strength constant within a distance of about 1 meter. The reflected PPF density n^^^ will 
be received by a microwave receiver as a signal of GWs. As a typical pattern, some of the 
generated PPFs will travel around the maser beam [16], forming a circular flux, which is shown 
in Fig. 2. 

Although the maser beam is set to propagate along the positive 2;-direction overwhelmingly, 
there is always a leakage flux of photons along the radial direction (normal to the z-direction). 
These leaking photons will mix up with the perturbed photons and form a noise for detection. 
Since we have chosen to detect the PPFs in the x-direction, we should try to eliminate the 
BPFs in the x-direction. This can be achieved by using a linearly polarized maser beam. Maser 
beams in the GHz band have been generated under laboratory conditions [20, 21]. In our setup 
of the detector, we make use of a maser beam with transverse electric mode (TEMqo), whose 
strength has a Gaussian-type distribution [35]: 

V 1 + (^//) '^ J ) 

where ■i/'o is the amplitude on the plane z = 0, r = ^x^ + y^ is the radial distance, W = 
Wq[1 + (2;//)^]^/^ with Wq being the radius of the beam on the plane z = 0, f = W^ke/^, 
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Figure 1: The sketch map of the setup of detection. The maser beam is propagating along 
the positive z-direction, and the static magnetic field points to the positive ^/-direction as 
represented by the furcations. The envelope of the maser beam is sketched as the hyperbola. 
The fractal membrane is placed in the y — z plane near the maser beam, facing to the right. The 
portion of the PPFs along the negative x-direction is reflected by the fractal membrane. The 
reflected PPF density n^^-* goes along the positive x-direction and keep its strength constant 
within 1 meter. The microwave receiver is placed on the right, facing to the left, and catching 
the outgoing perturbed photons as the signal. 




Figure 2: A sketch map of the circular perturbed photon flux. The inner circle stands for the 
maser beam. The dark areas stand for the perturbed photon flux in the circular direction, 
where the thicker region means a larger flux. The membrane is placed near the Gaussian beam 
and the flux is reflected back into the positive x-direction. 



ke = 27r/Ae is the wavenumber, Ue is the angular frequency, R = z + P j z is the curvature 
radius of the wave fronts of the beam at z, and b is an arbitrary phase factor. Furthermore, 
the maser beam is hnearly polarized along the x-direction, namely, the electric field in the 
maser beam is given by 



Ef(x,t) = ^(x,t), Ef = ^f = 0, 



(2) 



where the tilde "~" and the superscript "(0)" stand for the time-dependent and the background 
EM fields, respectively. Since the maser beam emitted from the emitter is in the region 2; > 0. 
the static magnetic field is chosen to localize in the region 2; > 0, 



m = I ^f 



{0< z< I), 

{z <0 or z > I), 



(3) 



where the caret "A" denotes the static magnetic field, and / is the dimension of the static 
magnetic field in the 2;-direction. 

In absence of GWs, the components of the average BPF density (m^'^s^^) are given as 
follows. 
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where "()" means the average over a time scale much longer than l/z/g- Note that, since 
the maser has been chosen to be linearly polarized in x-direction, the x-component of the 
BPF density, n}^' as a noise, is vanishing. This feature is an advantage over that using an 
unpolarized beam [19]. Of course, in actual situation, the maser beam can not be polarized 
completely. Then, Eq.(4) is not valid strictly, and it always exits the residual BPF density 
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fractal membrane only reflects the PPF not the BPF [19]. However, after reflected by the 
fractal membrane, ni^^ will not decay within 1 meter. Then, at a large radial distance r from 
the beam, n)^' can be negligible compared with n^J-'. We will discuss about this problem in 
more details in Sec III. 

III. PPFs GENERATED BY GWs ALONG VARIOUS DIRECTIONS 



A. GWs along positive ^-direction 



The situation we are interested in is when GWs are present. Consider a beam of circularly 
polarized monochromatic plane GWs propagating along the positive 2;-axis. The metric can 
be written as 

where r]^,^ is Minkowsky metric, and h^^ stands for GWs with \h^jy\ ^ 1. In TT (transverse- 
traceless) gauge, h^^ has only two independent components: hn = — /i22 = h^, and hi2 = 
h2i = h^, where 

/i® = A®exp [i{kgZ - ujgt)], 

h(x, = iAtg) exp [i{kgZ - cOgt)]. (8) 

In a curved spacetime, the Maxwell's equations in vacuum are [36, 37] 
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where Fy,y is the EM field tensor, g = — det((7^jy), and the comma means the ordinary derivative. 
Since the EM field will be perturbed by GWs, we decompose the total EM tensor into two 
parts: 

where FJ^^ represents the background fields, and FJ^^) the perturbed fields caused by GWs. 
Explicitly, 
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Since \hnu\ <^ 1, FJ^^ is also small and will be evaluated up to the first order of \huu\- As 



will be seen, each component F^^ receives two parts of contributions: one comes from the 
interaction between the static magnetic field and the GWs, oc \h^^\B^^\ the other comes from 
the interaction between the maser beam and the GWs, oc \h^^\B'^'^^ [38, 39, 15]. In our designing 
of the detection, the static magnetic field is chosen to be so large that B^^^ /IB'^'^^ ~ 10^^ [19]. 
Therefore, in F^]} we only keep the contribution oc \h^u\B^^\ The maser beam just provides 
the resonance condition to generate the PPFs, i.e., the detector only respond to the GWs with 
the same frequency as the maser beam. By solving the Maxwell's equations (9), one obtains 
the perturbed EM fields and the PPFs. The detailed calculation of n'^^ is given in Appendix. 
The resulting expressions of the PPFs density n'^^ in all three regions are given by 



Region I (-2 < 0). 

n« = 0; (13) 

Region II (0 < ^ < /), 

(1) ^^^S°Vol/ i kgz 
n^ ' = < Rin O 

Region III (z>l), 

2^„hu,.ii + (z/fy-v'^\2(z + fyz) 

where /xq is the permeability in vacuum, and 

<^ = ^-arctan(^). (16) 

The phase 6 = 7r/2 has been taken for concreteness. As Eqs.(14) and (15) show, n^p is only 
produced by the GWs of (^-polarization mode, and is proportional to the static magnetic 
field B^^^ and the maximal amplitude ^q of the maser beam. Since |n^^^| contains a decaying 
factor e~^ /^ , it decreases radially for larger r. To visualize the dependence of n^^^ on spatial 
variables, we plot it as a function of {y, z) on the plane x = 0.05 m in Fig. 3, and as a function 
of (x, y) on the plane z = 0.4 m in Fig. 4. 

B. GWs along some other directions 

The above is for the incident GWs along the positive 2;-axis. In this section, we give the 
results of incident GWs propagating along other directions, while the setup of the detector is 
the same as in Section II. 

(1) The incident GWs propagating along the negative z-direction. The relevant results are 
the following: 

Region II (0 < z < /), 
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Figure 3: n^^^ as a function of (y, z) on the plane x = 0.05 m. The parameters are taken as 
A^ ~ 10-3°, ipo = 1.8 X 10^ Vm-\ Wq = 0.05 m, S(°) = 3T, and / = 0.4 m for demonstration. 
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Figure 4: n^^-* as a function of (x, ?/) on the plane z = 0.4 m. The parameters are taken the 
same as in Fig. 3. 



Region III {z > I), 

n« = 0. (18) 

Since the maser emitter lies in Region I {z < 0), so it is not interesting for detection. 

(2) The incident GWs propagating along the positive x-direction. The static magnetic field 
B^^' is taken to be localized in the region —li < x < ^2 in the x-direction. One obtains 

^^ - 2^onu;.[l + {z/m^- \ 2{z + P/z) '^" ^^'^' "^ + ^^ 
X + li 1 

^ Wi[l + {z/ff] ''''^ '^^^"^ - a:) + $] + 2{z + p/z) ""^^ ^^^^^ ""^^ ^^'^ ^ ^^ 

^ k,wi[i\{z/m ^^" (^^") ^"^ ^^^' + ^) } ^"p ['w^ ('') 

for the region — /i < x <l2- 

(3) The incident GWs propagating along the negative x-direction. With the static magnetic 
field as in (2), one has 

n« - ^^BfH,y f h{x - h) 

^k^w![iTWm '" ^^^"^^ '°' ^^^' + ^^} "^ (~^) ' ^'°^ 

which is similar to Eq.(19). 

(4) The incident GWs propagating along the positive or negative y-direction. One finds 
that -EW = 0, leading to 

4^) = 0. (21) 

The above results show that, for the given setup, the detector responses differently to the 
incident GWs from different directions. In the following it will be seen that the detector 
responses most effectively to the GWs in the 2;-direction. In general, RGWs is of stochastic 
nature and come from various directions. In this case, a reduction factor will be introduced as 
shown later. 

C. Numerical calculations for PPFs 

In order to examine the dependence of perturbed photons on the directions in which GWs 
propagate, let us estimate numerically the perturbed photons per unit time received by the 
microwave receiver for the above cases of the incident GWs. For concreteness, we adopt the 
following parameters of the detector that can be realized in the laboratory: 



Table 1: The number of photons per second passing through As, A^^^-*, for GWs from various 
directions. 
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1) P = 2kW, the power of the maser beam, corresponding to ipQ = 1.8 x 10^ Vm ^ for the 
spot radius Wq = 0.05 m [20, 21]. 

2) B^^^ = 3Tesla, the strength of the background static magnetic field [40]. 

3) / = 0.4 m, the width of the static magnetic field in z-direction. 

4) I' = li + 12 = 0.4 m, the width of the static magnetic field in a;-direction. 

5) z/g — 4.5 GHz, the frequency of the maser beam in the microwave band [20, 21]. 

The number of perturbed photons in the x-direction per second passing through a surface 
As on the plane x = 0.05m is given by: 

Ni'^= fn^^^U=omdydz. (22) 



Here the integrand is taken to be the negative portion of n^^-* < 0, which is reflected by the 
membrane back to the positive x-direction. For comparison, we choose As ~ 8 x 10~^m^ 
(0 < y < 0.2 m, < 2; < 0.4 m) to receive more photons with a limited size. For concreteness, 
A^ ~ 10^^'^ is taken. The resulting A^^^^ is shown in Table 1. We see that the magnitude of 
iV|^) generated by the incident GWs along the positive 2;-direction and the positive x-direction 
has the same order, which is larger than that for other cases. Note, for the case of GWs along 
y-direction, N^^^ is nearly vanishing. 

Using the parameters given above, if we choose the maser beam to have a polarization 
degree is ~ 98%, corresponding to a ratio of unpolarized/polarized electric field components 
E^^^ /E^^ ~ 0.1, our computation shows that the ratio of the number of background/perturbed 
photons over the area As per second will be 

Ni'^/Ni'^ ~ 10-^ (23) 

at a; ~ 1 m . So, if the microwave receiver is put at ~ 1 meter away from the fractal membrane, 
the influence of the background photon flux will be effectively negligible. A higher polarization 
of the maser beam is always wanted to suppress the BPF. 

Remember that the phase 6 = 7r/2 has been taken in the above for simplicity. However, 
in general, A^^^^ would depend on the phase factor 5. Fig. 5 (a),(b),(c) and (d) give A^^^^ as a 
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Figure 5: The small fluctuation of N^'^ with variable 5. (a),(b),(c) and (d) give the cases of 
the incident GWs along +2;, — z, +x and — x, respectively. 

function of 5 for the incident GWs along +2;, — z, +x and — x, respectively. Fig. 5 shows that 
the changes of A^^^-* with 5 is small, and the error for various values of 5 is less than 4%. Thus, 
in the following, we assume that A^^^-* is independent of the phase factor 5. 

IV. THE DETECTION FOR RGWs 

In this section, we estimate the sensitivity of the detector, and analyze the feasibility of 
detecting the RGWs using this method. 

A. Implementation of the experiment 

One certainly expects to have a number of problems in the actual implementation of such 
kind of detection. There could be various sorts of noises for this detector, such as thermal noise, 
external EM noise, seismic noise, and shot noise in the maser beam, etc. The seismic noise 
is one of major obstacles for the common laser interferometer detectors. However, this kind 
of noise usually has a frequency much lower than GHz band and will not generate additional 
perturbed photons, so it will not affect our detection essentially. Still, an isolating system 
may be employed. For instance, the detector system may be put on a suspended framework 
to absorb seismic vibrations. As for the shot noise in the maser beam, it can be suppressed by 
stabilizing the frequency and the amplitude of the maser beam. 

Among all these kinds of noises, the external EM noise would be a real problem for our 
detection. For example, the CMB at T ~ 2.7 K yields a photon flux density n ~ 10^ cm^^s^^ 
around the frequency v ~ 4.5 GHz with a width Av ~ 10 KHz (the frequency width of 



the maser) in any direction. The energy flux of CMB photons would overwhelm that of the 
perturbed photons, since the average PPF density fS^^ ~ 0.6 cm~'^s~^ as from Table 1 where 
we have assumed A^ ~ 10~^°. To solve the problem, we propose to employ a Faraday cage, 
which shields the detecting device from the external EM noise. The outer shell of the cage 
may be made of some conducting metal so that the external EM waves will not entering the 
cage. The inner surface of the cage should be made of some kind of material that effectively 
absorbs noise photons within the cage. So the CMB photons and other EM noise inside the 
cage can be eliminated. In designing the Faraday cage, we should excavate a little hole on the 
cage and the hole should be sealed by a one-way membrane with a total transmittance around 
4.5 GHz, so that the maser beam can pass through it while the photons can not enter the cage. 
Furthermore, to eliminate thermal photons emitted from the detector system and from the 
inner layer of the cage, one need reduce the temperature of the system. Therefore, a cryogenic 
technique should be applied so that the detector operates in a low temperature environment. 
Moreover, a vacuum environment of the system will help for the detection. 

B. Sensitivity of detector 

For a preliminary analysis, we will focus on the thermal noise in our detection system and 
estimate its sensitivity limited by thermal noise. The signal power is given by 

^ = r/|iV«|^a;„ (24) 

where r] is the reflectance of the fractal membranes, ranging from to 1. With the help of 

Eq.(22), iVW is given by 

HohuJe 

The input thermal noise (the thermal noise coming into the input part of the receiving system) 
can be estimated as 

iVin = keTB, (26) 

where ks is Boltzmann's constant, T is the temperature of the thermal noise, and B is the 
bandwidth of the detector in Hz, which can be estimated as -B ~ Ve/Q-, where Q is its quality 
factor. There are additional thermal noises within the receiving system, thus the minimal 
signal power should satisfy [42] 

5mi„ = M/oiVi„, (27) 

where M > 1 is the minimal output signal-to-noise, and /q > 1 is the noise coefficient of the 
microwave receiver, defined as the ratio of the input signal-to-noise to the output signal-to- 
noise. Then using Eqs.(25)-(27) and letting M = 1 yield the minimal detectable dimensionless 
amplitude, 

ft,,,,. ^ <'Jt^ X 10". (28) 



where /i^in = ^^min- Taking ^o = 1.8 x 10^ Vm-^ and 5j°) = 3 Tesla, Q = 10^ T ~ ImK 
[43], /o ~ 2, ?7 ~ 99.99% [22, 23, 24], one obtains the sensitivity 

/iinin ~ 4 X 10~^°. (29) 

As said earher, RGWs come from all directions and form a stochastic background, therefore, 
in evaluating the sensitivity of the detector, a reduction factor F should be introduced [41]. We 
can estimate its magnitude as follows. Firstly, we consider the case that the detector responses 
only to the incident GWs along the positive 2;-direction. By Eqs.(42) and (57) in Appendix, 

n W oc Ef oc dh^ {z, t) jdz oc kji^ {z,t). (30) 

For a beam of incident GWs with a wave vector k, one needs to project it along 2;-direction, 
so that its component n]^' oc kz = kcosO, where 9 is the angle between k and the positive 
z-direction. Then in this special case the reduction factor F will be estimated as 

F = ^ f\os9smede T" d(p=\. (31) 
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However, this estimate is not complete. In fact, while the detector does not response to 
the incident RGWs in the y-direction, it responses to the incident RGWs in both the x- 
and 2;-directions with the same order of magnitude, as shown in Table 1. For an qualitative 
estimation, we can take the response of the detector to any incident GWs perpendicular to the 
^/-direction to be the same. Consider an arbitrary beam of GWs whose wave vector k forms 
an angle 6 with the y-direction. Projecting the wave vector k on the x — z plane gives rise to 
a factor ~ k sin 6. Then, taking the average of sin 6 over the solid angle Air yields 

F=— r sin^ 9d9 /^" rf0 = -. (32) 

Att Jo Jo 4 

Therefore, we expect that the actual reduction factor F for our detector would be between 1/4 
and 7r/4. Multiplying Eq.(25) by F, the sensitivity of the detector given by Eq.(29) should be 
modified as 

/i„,i„ ~ (5.1 X 10"^° ~ 1.6 X 10"^^). (33) 

C. Detecting RGWs 

What about the detection target, say, the the RGWs in the present accelerating universe 
[27, 28, 32] around the frequency Ug ~ 4.5 GHz? Now we calculate the root-mean-square 
(r.m.s.) amplitude of RGWs. In the high frequency limit, the RGWs can be considered 
approximately as the superposition of plane waves in Eq.(8). By its nature, RGWs constitute 
a stochastic background, and the mean value of the field hij is zero at every instance of time 
and at every spatial point: (0|/ijj(x, r)|0) = 0. But the variance is not zero [25, 27, 28], 
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Figure 6: The spectrum of RGWs for the cosmological model with the tensor-scalar ratio 
r = 0.22, the dark energy Q\ = 0.75, and the inflation parameter (3 = —1.9. The spectrum in 
GHz band depends sensitively on the reheating parameter Ps [29]. 

where /i(z/, r) is the spectrum of the RGWs, and r is the conformal time. Fig. 6 gives the 
spectrum /i(z/, th) of the RGWs at the present time th for the cosmological model with the 
tensor-scalar ratio r = 0.22, the dark energy Q^ = 0.75, the inflation parameter jS = —1.9, 
and the reheating parameter Ps [29]. The quantity hiy^Tu) is related to the spectral energy 
density VLgiy) often used in literatures [25, 29, 41], 
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where z/j;/ = ifo ~ 2 x 10^^^ Hz is the Hubble frequency. 

Due to the resonance condition, the detector only responses to a very narrow frequency 
band Az/ ~ Vg/Q around the central frequency Vg, where Q is the quality factor of the maser 
beam. Thus, only the modes of frequencies Vg ~ 4.5 GHz are selected among the incident 
RGWs. The integration in Eq.(34) is then evaluated as 
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So the r.m.s. amplitude of the RGWs in the band is 
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(37) 



where \/2 accounts for the assumption that the ®- and ©-polarization modes give equal con- 
tribution. Reading from the known spectrum in Fig. 6 gives 



h{vg, Th) ~ 8 X 10 



-32 



(3^ 



at z/ ~ 4.5 GHz for the reheating model I3s = 0.5. The corresponding r.m.s. amphtude is 

/irms ^ 1-8 X 10"^^ (39) 

Comparing Eq.(33) with Eq.(39), one can see that there is approximately a gap of 4 ~ 5 orders 
of magnitude between the sensitivity of the sample detector and the r.m.s. amplitude of RGWs 
in the accelerating universe. 

As for the possible thermal background of gravitational waves [33, 34], it could have been 
generated at the very early stage at an energy scale ~ 10^^ Gev, also described by the Planck 
spectrum like CMB photons. If there is no inflationary process, the graviton gas would be at 
T ~ 1 K, corresponding to typical frequencies in ~ 100 GHz. The amplitude of the spectrum 
of this thermal GWs would be about /i(z/) ~ 10~^^ around the frequency of 4.5 GHz. However, 
if the inflationary expansion has occurred by some 60 e-folding, as is supported by WMAP 
data of CMB anisotropies [44] and others, the thermal GWs would be drastically diluted and 
its temperature would be reduced to T ~ 10~^^ K, totally negligible. 

D. Possible improvements 

Although the sensitivity is still 4 ~ 5 orders short to detect the RGWs, the detector has 
a large room for improve in several ways. Firstly, note that /imin oc 1/Q by Eq.(28), while 
/irms oc 1/v^ by Eq.(37), so the ratio Kms/Kmn oc -/Q. A larger quality factor Q of maser 
beam will enhance the possibility for detection. This would require a highly monochromatic 
maser beam. For instance, if Q can be increased from 10^ to 10^, the ratio would be enhanced 
by ~ 100 times and the gap will be reduced by 2 orders. At present, for conventional lasers in 
the optical frequency band, a quality factor Q ~ 10^^ has been achieved [45], and for hydrogen 
maser, the quality factor Q has been reached up to ~ 10^ [46, 47, 48]. Secondly, as can be 
seen from Eq.(28), the sensitivity depends strongly on the temperature T of detector. If it 
is reduced down to T ~ 50/iK [49], the sensitivity will be improved by a factor ~ 20 and 
the gap will be suppressed to be about 1 order. Thirdly, increasing the strength of the static 
magnetic field B^^ and the power of the maser will also improve the sensitivity of the detector. 
Apart from the above possible improvements, enlarging the interaction dimension between the 
GWs and the static magnetic field will also improve the detection [19]. Putting these possible 
improvements together, an actual detection will be realistic. 

V. CONCLUSION AND DISCUSSION 

We have extensively studied the maser beam method for detection of GWs ~ 4.5 GHz. The 
experimental setup consists of a maser beam, a strong static magnetic field, a reflecting fractal 
membrane and a microwave receiver. Moreover, a Faraday cage should be used to prevent 
the detector from external EM noises. And, to reduce thermal noise, the detector should be 
place in a low temperature environment. The maser beam is chosen to be linearly polarized. 



so that the BPF in the detecting direction can be suppressed effectively, and the PPF can be 
detected as a signal of GWs. We have obtained the analytical expressions for the PPF density 
n^^-* generated by the incident GWs from various directions. 

To examine the feasibility of the detection, we have estimated the sensitivity of the sample 
detector limited by thermal noise and have confronted it with the RGWs in the accelerating 
universe as a scientific object. In our preliminary analysis, we found that there was still a 
gap of about 4 ~ 5 orders between the sensitivity of the detector and the r.m.s. amplitude of 
the RGWs. However, we have a lot of ways in improving the sensitivity, such as lowering the 
temperature, increasing the quality factor and the power of maser beam, and enlarging the 
strength and dimension of the static magnetic field. These improvements will remove the gap, 
making the method applicable for detecting high frequency RGWs. 

However, our analysis on the detection for the RGWs are still tentative, and the conclusions 
arrived are also preliminary. In particular, the odds is that the detecting PPF as the signal is 
very small and confronts a number of possible sources of EM noise. A systematical analysis 
on effectively suppressing these noises is then needed to give a more reliable sensitivity. 

Overall, the maser beam method in GHz band or higher is feasible. As a new method 
to detect GWs, it is complementary to the laser interferometer method working in the low 
frequency range (10^"^ ~ 10^) Hz. Moreover, from the point of view of experimental construc- 
tions, the building of this detection is much less expensive than ordinary interferometer laser 
methods. Therefore, under these considerations, the GW detection scheme is certainly worthy 
of further studies and is expected to be implemented in laboratory someday. 
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APPENDIX: CALCULATIONS OF n^ FOR GWs ALONG DIRECTION OF 

MASER BEAM 

In this appendix, we present the calculations of the PPFs density n^p produced by the 
GWs along the positive ^-direction. By Eq. (3), the space is divided into three regions: I 
{z < 0), II {0 < z < I), and HI {z > I). Firstly, we focus on the region II, where the static 
magnetic field B^^^ ^ 0. By Eq.(8), the expressions of GWs only have two variables {z,t), so 
will be the perturbed EM fields accordingly. Plugging Eqs. (7) - (12) into Eqs. (9) and (10), 
and keeping only up to the linear terms of hn„, then after lengthy but easy calculations, one 
obtains 



1 ^ 
c 



^E« + 5(;)=5Wv. (40) 



TZ> 

y,z 
and 



^i:i + B^^t = (41) 

^Xl-Bi]l = Bi'^h^,. (42) 

m - m = 0, (43) 



Ei]l = Ei]l = 0, 5«=5«=0. (44) 

Using Eq.(8), one solves Eqs.(40)-(43) in Region II [16]: 

E« = -A^B^^hgZe'^'"''-'^'*'' + 62e*(^'«^"^«*) + c2e^('=«^+'^«*), (45) 

E^i) = iA^5(°)A;^2e^('^«^-'^«*) + i64e*(^«^-'^«*) + ic4e'^''^'+'^^''> . (46) 

From Eq.(44), one obtains a physical solution, 

Ei'^ = B^p = 0, (47) 

which is valid in all the three regions. The constants, bi, Ci, ... 64, C4, in Eqs.(45) and (46) are 
to be determined by the physical requirements and boundary conditions in the following. 

Any physical measurement by an observer in curved spacetime should be carried out in a 
local inertial frame, i.e., the observable quantities are the projections of the physical quantities 
on to the four orthonormal bases e^, e^, e^, e^ carried by the observer. Therefore, the observable 
EM fields are 

Fa0 = F^.e^^^'f,. (48) 

For an observer at rest with respect to the static magnetic field, one can choose 



e[( = (1,0,0,0) 
1 
2 



^ = (0,l--/ie,0,0). 



e^ = (0,-/i«,l + i/ie,0), 

e^ = (0,0,0,l). (49) 

Suppose that, for simplicity, there is no perturbed EM waves propagating in the negative z 
direction in Region I and Region III [16]. From Eqs.(12), (48) and (49), and by the boundary 
conditions that the real parts of the perturbed fields F^J are continuous at the interfaces, the 
observable perturbed EM fields in the three regions are given by: 
Region I (5(0) = 0): 

Ff = Ef' = B^ = Bf' = 0; (50) 



Region II (5^ = 5^): 



eP = ^Ae5(°)cA;gZe*(^«^-^«*), 



2 ^ y 9 
-A r(°)/^ ' 



Bf = -Ae^f A;,^e*('^«^-"''*), (51) 



5(1) — _A R(0)i, 7p*(fc9Z-<^9t) 

+^A^M^e'^''^'-''^'^ + ^A^M^^e'^^^'+'^^'l (52) 

Note that the expressions of -E^ and B^ in Eq.(51) are different from those given by Eq.(43) 
inRef.[16]. 



Region III (fi(o) = 0): 



Bf = ^AeSf fc,/e*(^«^--«*), (53) 

5^ ) = ^A^5f A;g/e*(^«^-'^«*), (54) 



where / satisfies 



/ = n\g (n is an integer). (55) 

It is straightforward to obtain from Eq.(47) 

Ef = 5« = 0, (56) 

also valid in all the three regions. From the perturbed EM fields given above, it is straight 
forward to obtain the PPFs density: 

n« = ^(4'^^f )^.-e, (57) 

where the subindex "z/^ = z/e" indicates the resonance condition, under which the time average 
will be non- vanishing. 
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